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jl.  Introduction 


1  In  this  paper  we  present  results  related  to  the 
smoothing  problem  and  related  generalized  Riccati 
equations  for  the  two-point  boundary  value  descriptor 
system  (TPBVDS) 

Ex(k+1)  =  Ax(k)  +  Bu(k) 

Vj.xfO)  +  Vfx(N)  =  v 

j  y(k)  =  Cx(k) 

where  E,  A,  and  are  possibly  singular  nxn 

matrices,  and  B  and  C  are  nxm  and  pxn  matrices 
respectively. 


(1) 

(2) 

(3) 


]  II.  System '  Theory  for  TFBVDSs 


I  In  [1-2]  we  develop  a  basic  theory  for  (l)-(3). 
]Many  of  the  aspects  of  this  theory  have  a  similar 
flavor  to  that  in  [4-5],  although  the  possible 
singularity  of  E  and  A  creates  some  significant 
differences.  As  discussed  in  [1,2],  when  (l)-(2)  is 
iwell-posed,  we  can  assume  that  it  is  in  standard  form, 
ji.e.  for  some  constants  a  and 

aE+ph  =  T  (4) 

land 

V.EN+VfAN=I  .  (5) 

!  As  in  [4-5] ,  x(k)  can  be  decomposed  into  an 
'.outward  process  zq  and  an  inward  process  .  The 

outward  process  zq  is  defined  as 

zo(k,t)  =  Et-kx(t)  -  Ai-kx(k)  k<4.  (6) 

By  eliminating  x’s  in  (6),  we  find  that  z0(k, l)  is 
only  a  function  of  the  inputs  inside  the  interval 
[k,  t] .  Also  note  that  zq  does  not  depend  in  any  way  on 

the  boundary  matrices  and  V^.  The  expression  for 

the  inward  process  z^  is  in  general  very  complex, 

although  in  the  so-called  stationary  case  there  is  a 
simple  expression  for  z^  [1]. 

The  system  (l)-(2)  is  strongly  reachable  on  [k.t] 
jif  the  map  from  {u(m):  me[k,t-l]}  to  zQ(k,t)  is  onto. 

System  (2.1)  is  called  strongly  reachable  if  it  is 
i reachable  on  some  [k, 1] . 

iTheorem  1: 

!  The  following  statements  are  equivalent 

a)  System  (l)-(2)  is  strongly  reachable. 

b)  The  strong  reachability  matrix 

R  =  |An-1B !EAn_^B | • • • |En_1B j  (7) 

has  full  rank. 

c)  The  matrix  [sE-tA!B]  has  full  rank  for  all 

!(s,t)/(o,o). 

Id)  The  state  x(i)  where  ie[n,N-n]  can  be  made 
arbitrary  by  proper  choice  of  the  inputs 
u(j):  je[i-n, i+n-1]  with  all  other  inputs  and  the 
boundary  value  v  set  to  zero,  and  for  all  pair  of 
^matrices  Vj  and  Vf  in  standard  form. 
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i  The  system  ( 1 ) — (3)  is  strong ly  observable  on 
:[k,t)  if  the  map  z^k.i)-*  (y(m):  me[k,l]}  is  one  to 

ione.  System  (l)-(3)  is  called  strongly  observable  if 
it  is  observable  on  some  [k.t]. 


Theorem  2= 

The  following  statements  are  equivalent 

a)  System  (l)-(3)  is  strongly  observable. 

b)  The  strong  observability  matrix 

rCAn“1 
CEAn~2 


[CE11  1 
has  full  rank. 

; c)  The  matrix 


(8) 


sE-tA 

C 


has  full  rank  for  all 


,(s,t)*(0.0). 

d)  For  all  matrices  and  Vj.  in  standard  form,  the 


state  x(i)  where  ie[n,N-n]  can  be  uniquely  determined 
from  the  outputs  y(j):  je[i-n. i+n-1). 

It  is  also  possible  to  define  notions  of  weak 
reachability  and  observability  which  explicitly 
involve  the  boundary  matrices  and  and  to  develop 


la  theory  of  minimal  realizations  [1-2],  In  addition, 
in  [1]  we  develop  methods  for  the  recursive  solution 
of  (l)  and  develop  several  notions  of  stability  for 
TPBVDSs . 

III.  The  Optimal  Smoother 

Consider  the  system  (l)-(2)  together  with  the 
moise-corrupted  observations 

!  y(k)  =  Cx(k)  +  r(k)  k=l . N-l  (9) 

yb  =  ffix(°)  +  Vfx(N)  +rb  .  (10) 

Here  r(k),  r^.  u(k),  and  v  are  mutually  independent, 
;rb  is  a  zero  mean,  Gaussian  random  vector  with 


covariance  and  r(k)  is  a  zero  mean  white  Gaussian 


inoise  process  with  covariance  R. 

!  It  can  be  shown  [3]  that  the  smoothed  estimate 
ix(k)  satisfies  the  following  TPBVDS 


,[x(k+l)] 

lMk+l)J 


£$ 


x(k)l 

Mk)J 


f  o_ 

Lc'R 


y(k) 


j  .  k=l . N-l 

(11) 


where 


land  where  7. ,  7^  and  K  are  complicated  matrices. 


To  compute  the  estimate  we  can  use  any  of  the 
recursive  algorithms  developed  in  [1-2].  One  of  these 
is  the  so-called  two-filter  solution  in  which  the 
TPBVDS  dynamics  are  decoupled  into  forward  and 
backward  recursions,  followed  by  a  correction  to 
account  for  the  boundary  conditions.  A  necessary,  but 
not  sufficient,  condition  for  stability  of  a  TPBVDS  is 
that  it  is  forward-backward  stable,  i.e.  a  decoupling 
transformation  can  be  found  so  that  the  forward  and 
backward  recursions  are  both  stable. 


In  the  case  of  the  optimal  smoother,  it  is  shown 
in  [3]  that  if  the  following  generalized  Riccati 
equations 

6  =  A’ (E0-*E*  +  BQB')-1A  +  C’R-1C  (14) 

#  =  A(E’*-1E  +  C’R~^C)_^A'  +  BOB’  (15) 

have  positive  definite  solutions  'i  and  S  then  there 
exist  invertible  matrices  M  and  N  such  that 


MSN-1  = 

f1  -1°  -ll 

[0  A’S  XE»  J 

(16) 

r  ""1  .“"1  l 

IWN-1  = 

Tat  V*  ol 

L  o*  ij- 

(H) 

IMoreover,  the  eigenvalues  of  AT  E'4  and  A'S  E6 
iare  inside  or  on  the  unit  circle.  Equation  (3.5)  is 
j cal led  the  descriptor  Hamiltonian  equation  and  the 
above  decomposition  is  the  descriptor  Hamiltonian 
!  -1  -1 
' diagonal ization.  Of  course,  we  would  like  AT  E'4' 

^and  A*S  *E0  *  to  be  strictly  stable.  This  occurs  only 
when  the  descriptor  Hamiltonian  has  no  eigenmodes  on 
the  unit  circle  i.e.  it  is  forward-backward  stable. 


Theorgir-2: 

;  If  the  system  is  forward-backward  detectable  and 

stabilizable  (i.e.  the  modes  on  the  unit  circle  are 
strongly  reachable  and  strongly  observable)  then  the 
j  corresponding  smoother  is  forward-backward  stable. 

I  IV.  Generalized  Riccati  Equations 


In  this  section  we  study  the  generalized 
algebraic  Riccati  .equation. 


A(E>-1E  +  C’R-1C)-1A' 


+  BOB' . 


(18) 


Theorem  4: 

If  (E,A,B)  and  (C,E,A)  are  strongly  reachable  and 
observable  respectively  then  (18)  has  a  unique 
positive  definte  solution. 

The  approach  used  to  prove  this  theorem  is 
similar  to  that  in  [6]  for  the  standard  Riccati 
equation.  Details  will  be  presented  in  a  future 
paper.  Existence  proceeds  as  follows.  From  Theorem  3 
and  the  fact  that  eigenmodes  of  the  smoother  occur  in 


(19) 


'  E 

-BQB’' 

F 

A  O' 

F 

.  0 

-A' 

G 

““ 

-C'R~*C  -E’. 

G 

I  The  proof  then  proceeds  by  first  showing  that  F  is 
then  that  E'GF 


•  invertible, 
that 


*  +  C'R  >  0  and  finally 


E*— 1 


+  C’R 


+  BOB’); 


<#>  =  (A(E’GF 
satisfies  (18). 

To  prove  uniqueness,  let  <p±  and  <pz  be  two 
positive  definite  solutions  of  (18),  let  Af  = 
and 


(20) 


Tt  =  EV^E  +  C’R_1C  for  i=l,2.  (21) 

Some  algebra  then  yields 

A?  =  AT^EV^A-p^ET'V.  (22) 

—1  —1  -1—1  * 

But  AT^  E.'ipj  and  ETg  A'  are  strictly  stable,  (see 

!  [3]);  thus  A<f  =  O'. 
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